Solidity of viscous liquids. III. Alpha relaxation 



in 
o 
o 

(N 



o 

c3 



c 

o 
o 



> ■ 
r- ■ 

: 

\o . 

(N ■ 

co : 

o . 

-i— > ■ 
a ■ 



i 

-a 
c 

o 
o 



X 
S3 



Jeppe C. Dyre 
Department of Mathematics and Physics (IMFUFA), 
"Glass and time" - Danish National Research Foundation Centre for Viscous Liquid Dynamics, 
Roskilde University, Postbox 260, DK-^000 Roskilde, Denmark 
(Dated: February 2, 2008) 

It is suggested that the a; -1 / 2 high-frequency decay of the alpha loss in highly viscous liquids, 
which appears to be generic, is a manifestation of a negative long-time tail as typically encountered in 
stochastic dynamics. The proposed mechanism requires that the coherent diffusion constant is much 
larger than estimated from the alpha relaxation time. This is justified by reference to the solidity 
of viscous liquids in an argument which, by utilizing the irrelevance of momentum conservation at 
high viscosity and introducing a center of mass diffusion constant, implies that at high viscosity the 
coherent diffusion constant is much larger than the incoherent diffusion constant. 



I. INTRODUCTION 

Viscous liquids and the glass transition have attracted 
considerable interest for many years 0, 0, El IE El- 
The solidity of viscous liquids refers to their solidness, 
i.e., that a viscous liquid is similar to a disordered solids 
because molecular motion is mostly vibrational. This 
was the view, e.g., of Goldstein who in his now classic 
1969 paper expressed the belief that "when all is said 
and done, the existence of potential energy barriers large 
compared to thermal energy are intrinsic to the occur- 
rence of the glassy state, and dominate flow, at least 
at low temperatures" Q. This implies that solid-state 
elasticity arguments are relevant and useful for under- 
standing viscous liquid dynamics 0. Actually, the view 
that viscous liquids are like solids by mostly being in 
elastic equilibrium is obvious already from Kauzmann's 
1948 review, which referred to these large barrier transi- 
tions (flow events) as "jumps of molecular units of flow 
between different positions of equilibrium in the liquid's 
quasicrystalline lattice" In the same spirit, Mooney 
in 1957 wrote that "a liquid not only could be but per- 
haps should be treated as an elastic continuum with a 
stress relaxation mechanism" p|. 

Based on the solidity property a theory for the non- 
Arrhenius temperature dependence of the viscosity was 
proposed some time ago, according to which the acti- 
vation energy is proportional to the instantaneous shear 
modulus (the shoving model) 0. This was justified by 
arguing that during a barrier transition the surrounding 
liquid behaves like a solid subject to elastic deformation, 
and that most of the barrier energy is shear elastic energy 
located in the surroundings. In two subsequent papers (I 
and II, E3) we further discussed the solidity of viscous 
liquids and introduced a "solidity length," I, giving the 
length scale below which a viscous liquid is virtually in- 
distinguishable from a disordered solid - albeit one that 
flows. 

By proceeding inductively from experiment the present 
paper argues that solidity may explain generic features of 
the primary relaxation process in viscous liquids, the so- 
called alpha process. Alpha relaxation is the slowest and 
the dominant relaxation process. The alpha process is re- 



sponsible for the glass transition, the kinetic falling out of 
equilibrium taking place when the coolingrate dlnT/dt 
exceeds the alpha relaxation rate 0, 0,13) IE El • Alpha 
relaxation is observable in dielectric relaxation experi- 
ments and in measurements of the frequency-dependent 
shear and bulk moduli I3.|ll| and, e.g., of the frequency- 
dependent specific heat [13, UM ■ The loss peaks obtained 
by measuring these linear response properties are of sim- 
ilar shape, but usually not identical. The loss peak fre- 
quencies are also similar, but not identical, and their tem- 
perature dependencies are usually quite similar. In the 
following we shall not dwell on the dissimilarities. The 
focus is on the shear modulus, adopting the point of view 
that this is the primary response function while, e.g., the 
dielectric loss is determined by the coupling of molecular 
rotation to shear flow Q EE 111 111 El 

The characteristic features of alpha response functions 
x(w), as reported in the literature, are: 1) The loss peak 
frequency Lo m is non- Arrhenius with an activation energy 
which increases as temperature decreases; 2) on the low- 
frequency side of the loss peak the loss is virtually Debye 
like (oc u>); 3) on the high-frequency side there are signif- 
icant deviations from Debye behavior, here the loss de- 
cays following an approximate power law: oc tu~ n 
(n < 1). In many cases a wing is observed a few decades 
above the loss peak frequency where n changes to a lower 
value. 

The vast majority of experiments report values of n be- 
tween 0.3 and 0.7 [jjj- The exponent n is often tempera- 
ture dependent, implying a violation of time-temperature 
superposition. At the 1997 3rd International Discussion 
Meeting on Relaxations in Complex Systems it was sug- 
gested by Olsen that beta processes may play a role at 
much lower frequencies than gene rally expected, and that 
these could explain the wing [2(| • This was subsequently 
confirmed by Lunkenheimer and co-workers who, upon 
long-time annealing of propylene carbonate and glycerol, 
found that the wing indeed develops into a separate re- 
laxation process |2y. While the matter is still actively 
discussed, in the opinion of the Roskilde group the fol- 
lowing picture emerges |22] : Once the effects of beta pro- 
cesses are minimized by going to sufficiently low temper- 
atures (still in the equilibrium liquid phase), the generic 
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features of alpha relaxation come to light. Generic alpha 
relaxation obeys time-temperature superposition, and for 
the high-frequency asymptotic behavior the generic al- 
pha process is characterized by the universal exponent 
n = 1/2: 

Finally, generic alpha relaxation does not have univer- 
sal loss peak widths. This conclusion was reached from 
dielectric relaxation measurements on several molecular 
liquids |22|. but it applies also for the shear modulus 
G(uS) j23,|2J|, although shear mechanical measurements 
are considerably less accurate than dielectric measure- 
ments. 



II. THE BEL MODEL AND ITS FORMAL 
GENERALIZATION 

Curiously enough, the conjecture that n = 1/2 is a 
"fixed point," which equilibrium viscous liquids approach 
at sufficiently low temperatures, represents a step back- 
wards in time. Around 1970 it was noted that n is often 
close to 1/2, and there was a considerable interest in 
explaining this theoretically. Some models, which over 
time have been suggested to explain Eq. (JTJ or its equiv- 
alents, are: In 1948 Zener in his theory of "relaxation 
by thermal diffusion" for ultrasonic loss predicted that 
the internal friction at high frequencies varies as cj -1 / 2 
[25|. Glarum (in 1960), Doremus (1970), Phillips, Bar- 
low, and Lamb (1972), Bordewijk (1975), Kimmich and 
Voigt (1978), Wyllie (1979), and Condat (1989) all ar- 
rived at the w _1//2 -tail from essentially one-dimensional 
diffusion models [H El IS |H H iH E2 ■ This is also 
basically the mathematical mechanism in the model of 
Isakovich and Chaban, who in 1966 predicted a limit- 
ing high-frequency sound absorption per wavelength in 
highly viscous liquids varying as lj^ 1 / 2 by regarding the 
liquid as a micro-inhomogeneous medium |33j (see also 
Ref. 34]). In 1967 Barlow, Erginsav, and Lamb, in 
what is referred to as the BEL model, took an electri- 
cal engineer's approach and regarded the shear mechan- 
ical impedance of a viscous liquid as a parallel combina- 
tion of a Newtonian liquid and a Hookean solid, imply- 
ing an w _1 / 2 -tail for the high frequency shear modulus 
0T|3^,|33- Montrose and Litovitz in 1970 derived the 
w _1 / 2 -tail in a model invoking diffusion and relaxation 
of some unspecified order |37| . Majumdar, working in 
the time domain, arrived at the short time equivalent of 
the high-frequency tail in 1971 from the assumption that 
diffusive modes exist with amplitude oc k~ 2 |38j. Cunat 
in 1988 from a purely thermodynamic treatment arrived 
at a distribution of relaxation times on the logarithmic 
scale r proportional to -y/r for r — > 0, which implies the 
^-i/ 2 high-frequency tail [3|j. Sjogren in 1990 showed 
that Eq. JTJ arises in mode-coupling theory under certain 



conditions 0. Cichocki and Felderhof m 1994 de- 
rived the a; -1 / 2 high-frequency tail in a three-dimensional 
diffusion model with a radial potential barrier, which is 
mathematically similar to the Isakovich-Chaban model. 
- Most of these models involve some species obeying the 
diffusion equation, which is assumed to apply also at 
times shorter than the alpha relaxation time |42||. This 
nontrivial assumption is also central for the "long-time 
tail" mechanism discussed below. 

The BEL model focuses on the shear compliance 
J(u>) = 1/G(u>). In dimensionless units the BEL model 
has a definite prediction for J(u>), implying universal loss 
peak width. Because this prediction did not fit all experi- 
ments, the BEL expression was soon formally generalized 
by the same authors 01 to 

J(w) = a + &(iw)~ 1/a + c(iw) _1 . (2) 

The original BEL model corresponds to b 2 = 4ac. 
Clearly, Eq. J2J) implies Eq. (JIJ for the shear modulus loss 
G"(ui). In dimensionless units Eq. @, which we shall re- 
fer to as the "generalized BEL expression," has one free 
parameter and thus nonuniversal loss peak width. 



III. THE GENERALIZED BEL EXPRESSION AS 
A CONSEQUENCE OF A LONG-TIME TAIL 

Equation gives a good representation of generic 
alpha relaxation. We proceed to use the fluctuation- 
dissipation (FD) theorem to translate Eq. into a 
statement concerning a velocity autocorrelation function. 
According to the FD theorem, if Q is the x-direction dis- 
placement of the top of the liquid (area L 2 in the x-y- 
plane, height h in the z-direction) , the shear compliance 
is given by 

Here (AQ 2 (t)} is the equilibrium mean-square displace- 
ment of Q in time t, i.e., the zero shear stress thermal 
average of the square of AQ(t) = Q(t) — Q(0). 

Via Eq. the generalized BEL expression Eq. 
translates into 

(AQ 2 (t)) = A + Bt^ 2 + Ct. (4) 

In terms of the Q mean-square displacement the velocity 
autocorrelation function is given by 

<Q(0)QW> = \^(AQ 2 (t)). (5) 

The generalized BEL expression is thus mathematically 
equivalent to 
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(Q(O)Q(i)) cx - t^' 2 . (6) 

Equation ©, which is a simple way to summarize basic 
characteristics of generic alpha relaxation, immediately 
brings to mind the long-time tails discovered in 1967 by 
Alder and Wainwright in early molecular dynamics simu- 
lations of liquids . These authors showed that the ve- 
locity autocorrelation function for a single molecule does 
not decay exponentially to zero at long times as previ- 
ously expected, but as a power law with exponent —3/2: 

(x(0)x(t)) cx + t- 3/2 (t -> oo) . (7) 

The explanation is the following. According to statis- 
tical mechanics at any given time molecular velocities 
are uncorrelated. Consider a particular molecule at time 
t = 0. It has a certain momentum while the surround- 
ing molecules on average have zero momentum. As time 
passes the momentum of the molecule in question diffuses 
to the surrounding molecules (because the Navier-Stokcs 
equation is basically a transverse momentum diffusion 
equation). At time t the momentum has diffused the dis- 
tance ~ t 1 / 2 away, thus spread among ~ i 3 / 2 molecules. 
If the initial momentum of the molecule in question is 
roughly equally shared among these, on average every 
molecule has momentum ~ i~ 3 / 2 . This includes the par- 
ticular molecule in focus, thus explaining Eq. Q. - A 
nice classical review of long-time tails is that of Pomeau 
and Resibois from 1975 [4{|, while Kirkpatrick and co- 
workers recently gave an excellent review from a more 
general point of view |46| . 

Central to the above argument are momentum con- 
servation and the diffusion equation. If Eq. ijfJJl were to 
reflect a long-time tail, two questions come to mind: 1) 
What is the role of momentum conservation in viscous 
liquids? 2) How does one explain a negative long-time 
tail E3? 

The first question arises because viscous liquids close to 
the calorimetric glass transition have viscosities roughly 
10 14 times larger than that of liquids like ambient water. 
Since the kinematic viscosity is the transverse momentum 
diffusion constant, these liquids have extremely large mo- 
mentum diffusion constants. Consider a flow event taking 
the system from one minimum to another. If momentum 
conservation were relevant, a flow event could not result 
in movement of the center of mass. In reality the liquid 
is confined by container walls; these provide the external 
forces possibly needed to move the center of mass slightly 
to ensure that the molecules at the walls do not move. 
Thus just as momentum conservation plays no role for 
the description of point defect motion in crystals, mo- 
mentum conservation is irrelevant for the description of 
viscous liquid dynamics under laboratory conditions ^3 . 

Turning to the second question, it is shown in the next 
section that velocity autocorrelation functions are nega- 
tive in any stochastic dynamics, i.e., dynamics described 
by a master equation. 



What is the required conserved variable? Angular mo- 
mentum conservation is irrelevant for the same reason 
that momentum conservation is. Energy conservation 
could be relevant because heat conduction is rather slow, 
but this process is almost temperature independent and 
thus seemingly unable to give alpha relaxation times, 
which depend dramatically on temperature. The only 
alternative left is particle number conservation. 

Long-time tails arise from particle number conserva- 
tion because density fluctuations at long times decay fol- 
lowing the diffusion equation, leading to (where Ap is the 
deviation from average density) 

(Ap(r,0)Ap(r',t)> cx t~ 3 / 2 exp (- ^T") ' (8) 

From Eq. JSJl it follows that any generic variable F with 
zero average which is a [non-linear] function of density, 
decays following 

(F(O)F(t)) cx + t- 3 / 2 (t -» oo) . (9) 

The Appendix gives a simple example of how Eq. © 
implies Eq. ©. 

IV. NEGATIVE LONG-TIME TAILS IN 
STOCHASTIC DYNAMICS 

Viscous liquids are believed to be well described by 
stochastic dynamics ■ This could take the form 

of Langevin dynamics or, in the ultimate course-grained 
version, hopping between energy ("inherent dynamics" 
[5(j). We proceed to show that the velocity autocorrela- 
tion function in any stochastic dynamics can be written 
as the negative of a state autocorrelation function. If 
Pi is the probability for state i, the master equation is 
Pi = J2j A-ijPj- If t ne Green's function - the probability 
to go from state i to k in time t - is denoted by Gi^k(t), 
and -P C q,i is the equilibrium probability of state i, Eq. (J5J 
implies 

(Q(O)QW) = l^iQi-Qj^P^i^G^jit). (10) 

Using the fact that the Green's function solves the mas- 
ter equation {Gi^j = Ylk AjfcGi_+fc), detailed balance 
(P eq ,iGi^j(t) = P cq jGj^i(t)), and probability conser- 
vation (Y^i^-ij = 0), the right hand side of Eq. (|l(Jfl is 
reduced as follows: 

i,j,k 

= - y^(Q» — Qj) 2 A-jkP C q,kGk^i(t) 
i,j,k 
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= 2 ^ (Q * + Q i " 2 ftOj) A jfc p e q ,fcA li G fc ^ i (0 

i,j,k,l 

= - £ (QiA 4I )(QjA ifc )P eq , fc G fc _ > i(t). (11) 

i,j,k,l 

Thus in terms of the variable Bj* = ^ ■ QjAjk we have 



(Q(Q)Q(i)) = - <5(0)i?(t)> . (12) 

Since autocorrelation functions of any state function are 
always positive in stochastic dynamics, this proves that 
the velocity autocorrelation function is always negative. 
Moreover, as argued above, generic autocorrelation func- 
tions decay as i~ 3 / 2 for t — > oo. We thus get the required 
Eq. ©. 

V. CONDITION FOR LONG-TIME TAILS TO 
OPERATE AT RELATIVELY SHORT TIMES 



When asked what to expect if the coherent and inco- 
herent diffusion constants decouple, a majority of re- 
searchers opted for the opposite inequality. Equation 
(|15fl cannot be ruled out on experimental grounds, how- 
ever, because there seems to be no measurements of 
the coherent diffusion constant for liquids just above the 
calorimetric glass transition. 

VI. DIFFUSION CONSTANT DECOUPLING AS 
A CONSEQUENCE OF SOLIDITY 

In this section we prove that the required diffusion con- 
stant decoupling is consistent with the solidity of viscous 
liquids. More accurately, it is shown that if the coher- 
ent diffusion constant is frequency independent, solidity 
implies diffusion constant decoupling. We consider N 
molecules in volume V, however not with periodic bound- 
ary conditions. Density fluctuations are conveniently de- 
scribed in terms of the quantity pk defined [H^ by the 
following sum over all molecular positions Vj\ 



The negative t _3 / 2 -decay applies as t — > oo. The obvi- 
ous question is how a long-time tail could play any role at 
times shorter than the alpha relaxation time. The only 
possibility, it seems, is that density fluctuations, even at 
these relatively short times, decay follow the diffusion 
equation. Thus even for a few decades of time shorter 
than r a the coherent intermediate scattering function 
S(k, t) should obey 



S{k,t) 



cxp(-£>fc 2 i) . 



(13) 



If a is the average intermolecular distance, the coherent 
diffusion constant defines a characteristic time t c via t c = 
a 2 /D. It is natural to expect that Eq. (13) applies only 
at times much longer than t c . Thus the condition for the 
long-time tail mechanism to work above the alpha loss 
peak frequency is that t c -C r a , or 



(16) 



The static structure factor is given by S(k) = (/9kP-k) 
and the coherent intermediate scattering function by 
S(k,t) = (/9k(i)/0-k(O)) ,531]. In liquids and solids the 
long-wavelength limit of S(k) is a very small number 
equal to the ratio between system compressibility and 
ideal gas compressibility at same density and tempera- 
ture. Alternatively, the long-wavelength limit may be ex- 
pressed in terms of the number M of molecules in a large 
subvolume and its fluctuations: S(k -> 0) = (AM 2 )/{M) 
p3L l54j (which is one for an ideal gas by Poisson statis- 
tics, but much smaller for liquids and solids). This num- 
ber has two contributions, one from the phonon degrees 
of freedom and one from the inherent dynamics. The in- 
herent part of this small number |55| is thus even smaller. 
Henceforth S(k) refers to only the inherent part of the 
structure factor, so we clearly have 



— « D. (14) 

The single-particle [incoherent] diffusion constant D s 
is defined by D s = (Ax 2 (t)} jit (t — > 00). One naively 
expects that information about the structure is basically 
forgotten when - and only when - the molecules on aver- 
age have moved an intermolecular distance. This would 
imply that D s ~ a 2 /r a , which may well be the case for 
some liquids. In several liquids, however, it has been ob- 
served that D s is much larger than a 2 /r a 52]. In either 
case, Eq. (|14|) is clearly obeyed if the coherent and the 
incoherent diffusion constants decouple as follows: 

D a < D. (15) 



S{k^0) < 1. (17) 

We shall assume that inherent density fluctuations are 
described by a linear Langevin equation, 

Pk = -7(fc)Pk + aW- (18) 

The last term is the usual Gaussian white noise 
characterized by = £_k and (£k(t)£k'(*0) = 

27(fc)S'(fc)(5k+k',o<5(i — t'), which ensures that S(k) = 
(PkP-k)- We proceed to argue that in the k — > limit 
solidity implies that j(k) cx k 2 (as expected from parti- 
cle conservation) and D s <C D, where D is defined by 
j(k) = Dk 2 . In particular, Eq. (jT%|l implies Eq. (jT3|l for 
the intermediate scattering function. 
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It is assumed that I/ 1 / 3 is smaller than the solidity 
length. This implies that in the time between two flow 
events there is elastic equilibrium throughout the sample 
[To|. We may thus adopt inherent dynamics [Hcf. which 
regard flow events as instantaneous rearrangements tak- 
ing the system from one energy minimum (inherent state 
56]) to another. 

To prove that solidity implies ^(k) oc k 2 as k — > we 
first consider a single flow event. If Arj is the displace- 
ment of the j'th molecule, at small k the change in p^ is 
given by 

5 Pk = -^^e^ik-Ar,- (k-fO). (19) 



If ro marks the location of largest molecular displace- 
ments, at small k the exponentials may all be replaced 
by exp(ik • ro). Thus, if R is the sum of all particle co- 
ordinates we get 



oik-ro 



ik-AR (k -> 0) . 



(20) 



If momentum conservation were relevant, one would 
have AR=0, because the center of mass is fixed. A 
characteristic feature of viscous liquids, however, is the 
already mentioned irrelevance of momentum conserva- 
tion jig. Besides the coherent and single-particle dif- 
fusion constants there is consequently a third diffusion 
constant, which we shall refer to as the "center of mass 
diffusion constant," Dqm- In order to be well-defined in 
the N — ► 00 limit, writing R = (X, Y, Z) the center of 
mass diffusion constant should be defined by 



D c 



M 



(AX 2 (t)) 
2Nt 



(t 



(21) 



Here it is understood that the TV, V — > 00 limit is taken 
before t — > 00. If Dqm is frequency dependent, the above 
expression defines the dc limit. The high-frequency limit, 
-Dcm(oo), is determined by the mean-square displace- 
ment at short times: 



(|Ap k (t)| 2 ) = 2~f(k)S(k)t. Equating these two expres- 
sions leads to 



S(k -» 0) 



rfc 2 (Ax 2 ) 

2 7 (fc 0) 



(24) 



For the limit to be well defined and nonzero we must have 
7(fc) cx k 2 . 

Writing j(k) = Dk 2 , we proceed to show that solid- 
ity implies D s <C D. First note that, since Dcm(oo) = 
T{AX 2 }/2, Eq. (J2U becomes 



S(fc->0) = 



-Pcm(oo) 
D 



From the inequality (fTTj) we conclude that 



-Dcm(oo) < D . 



(25) 



(26) 



Consider a flow event taking the system from one min- 
imum to another. In the time between two flow events 
there is elastic equilibrium. The liquid is indistinguish- 
able from a disordered solids, and a flow event in- 
duces slight adjustments of the positions of molecules 
far from the relatively few molecules that move consider- 
ably 0, 1^3 • If the liquid were completely homogeneous 
and isotropic before and after the flow event, the total 
displacement AR would be zero. This is not realistic, 
however, so there is no reason to expect the total dis- 
placements to be zero. According to standard elasticity 
theory |10L 1.3 ll 1571 158| displacements far from the center 
of the flow event are small, varying with distance as r~ 2 . 
Their contribution to AR is insignificant, so we expect 



5> 2 ~ (5>, 



(27) 



This implies that the high-frequency limits of D s and 
Dcm are comparable: 



^cm(oo) = x (i -►()). (22) 

Over a short time span flow events are uncorrelated.Thus, 
if the number of transitions per molecule per unit time is 
denoted by T and (AX 2 ) is the average squared change 
of X in one flow event, the change in p^ in a short time 
t is given by 

(|Ap k (£)| 2 ) = NTt^k 2 (AX 2 ) = Ttk 2 (AX 2 ) (Jfc-fO). 

(23) 

On the other hand, the change in p k over a short time 
is determined by the noise term of Eq. (|18f) . leading to 



D s (w) ~ ^cm(oo). (28) 

Finally, we note that when diffusion constants are fre- 
quency dependent they always increase with frequency 
(see, e.g., Ref. |59(), and often quite a lot: 

D s < D s {oo). (29) 

To summarize, by representing the solidity of viscous 
liquids by inherent dynamics we have argued that 



— < D s < D a (oo) ~ £>cm(oc) < D . (30) 
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The strong inequality sign <C was used liberally through- 
out this section (the inequalities all hold with <C replaced 
by <, though): The first str ong inequality presumably 
applies only for some liquids |52|. while for others there 
may be near equality. The second strong inequality is 
likely - because this is what is generally observed for hop- 
ping in disordered solids |59j | - but not compelling. Only 
the third strong inequality is compelling, being a direct 
consequence of the small compressibility of liquids. Nev- 
ertheless, this is enough to establish the inequality (|14|l . 
once it is assumed that the coherent diffusion constant is 
frequency independent. 



VII. CONCLUSION 

A long-time-tail mechanism [60l | provides a natural ex- 
planation for the seemingly generic w~ 1//2 high-frequency 
decay of the alpha process. It has been argued that the 
required decoupling of coherent from incoherent diffusion 
follows from the solidity of viscous liquids. There are still 
problems to be addressed. One is that the derivation of 
the diffusion constant decoupling was based on the as- 
sumption that the sample is smaller than the solidity 
length. If the alpha relaxation time is one second, the 
solidity length is just a few thousand Angstroms [Io| . 
The question arises how to deal with bulk viscous liq- 
uids. Given the fact that Eq. I|18l) does not allow bulk 
volume relaxation on the alpha time scale, it is clear that 
something is missing j^. Nevertheless, we surmise that 
the long-time tail mechanism survives the necessary sup- 
plements to arrive at a full theory. Another problem 
is that the constants A and C of Eq. which here 

appear as simple integration constants, cannot be arbi- 
trary: In order for the generalized BEL expression to fit 
data these constants must be related to the constant B 
by AC ~ B 2 . A final challenge is to establish the con- 
sistency of the proposed scheme, where the coherent dif- 
fusion constant is assumed to be frequency independent 
whereas the incoherent diffusion constant most likely is 
strongly frequency dependent. 



APPENDIX: LONG-TIME TAILS DERIVED 
FROM THE DIFFUSION EQUATION 

We here show an example of how long-time tails typ- 
ically arise. By "long time" is meant times long enough 
that density fluctuations decay following the diffusion 
equation. As a simple case, consider a variable F which 
is a sum of pairwise contributions, F = J2ij 4>( r i ~ r j)- 
In terms of the density p(r) = J^i ^( r — r i) we have 



If the deviation from average density is denoted by Ap = 
p — (p) , F is given as 

F = Const. + / drdr' <j>(r - r') Ap(r)Ap(r') . (A.2) 
The ^-autocorrelation function is given by 

(F(Q)F(t)) = Const. + J dr 1 dr 2 dr 3 dr 4 ^(ra - r 2 )0(r 3 - r 4 ) 
x (Ap(ri,0)Ap(r 2 ,0)Ap(r 3 , t)Ap(r 4 , t)) 

(A.3) 

In the Gaussian approximation the average of a product 
of four variables by Wick's theorem is a sum of three 
terms, each a product of two pair correlation functions. 
One of the three terms is time-independent and we get 
[assuming </>(r) = </>(— r)] 

(F{0)F{t)) = Const. + 2 / dv 1 dv 2 dv 3 dv i (t){v 1 - r 2 )</>(r 3 - r 4 ) 
x (Ap(n, 0)Ap(r 3 , t))<Ap(r 2 , 0)Ap(r 4 , t)) 

(A.4) 

Letting f^oowe find (F) 2 = const, so, if AF = F—(F), 
Eq. 1A.4|) implies 



(AF(O)AF(i)) - 2 J dr 1 dr 2 dr 3 dr 4 0(r 1 - r 2 )0(r 3 - r 4 ) 

x (Ap(n , 0) Ap(r 3 , t)) (Ap(r 2 , 0) Ap(r 4 , t)) 

(A.5) 

A standard assumption in liquid theory is that the den- 
sity correlation function (Ap(r, 0) Ap(r' , t)) regarded as a 
function of r' and t obeys the diffusion equation at long 
times. Since the density correlation function goes to zero 
as |r — r'| — * oo and t — > oo, we have (where D is the 



diffusion constant) 



<Ap(r,0)Ap(r',i)) ex i" 3 / 2 exp (- ) . (A.(.) 



When this is substituted into Eq. I|A.5I) , one finds at long 
times, if 4> is reasonably short ranged, I = j dr<p(r), and 
V is the volume: 



F 



(AF(O)AF(t)) 



v i 2 r' i/2 . 



(A.7) 



drdr' • 



)p{v)p{v'). 



(A.l) 
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